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ABSTRACT 

Phases of the spherical harmonic analysis of full-sky cosmic microwave background 
(CMB) temperature data contain useful information complementary to the ubiquitous 
angular power spectrum. In this letter we present a new method of phase analysis on 
incomplete sky maps. It is based on Fourier phases of equal-latitude pixel rings of 
the map, which are related to the mean angle of the trigonometric moments from the 
full-sky phases. It has an advantage for probing regions of interest without tapping 
polluted Galactic plane area, and can localize non-Gaussian features and departure 
from statistical isotropy in the CMB. 
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1 INTRODUCTION 

The temperature anisotropy of the cosmic microwave back- 
ground (CMB) radiation contains a wealth of information 
about our Universe. Its statistical properties not only shed 
light on what kind of universe we are living in, but also 
lay the foundation for the significance and interpretation 
of the angular power spectrum. According to the gener- 
ally accepted cosmological model, namely the Cosmologi- 
cal Concordance Model, the primordial fluctuations in the 
early Universe constitute a Gaussian random field (GRF) 
(Bardeen et al. 1986; Bond & Efstathiou 1987). As the 
CMB is an observable imprint of the primordial fiuctuations, 
therefore, after the NASA WMAP data release (Bennett et 
al. 2003a; Bennett et al. 2003b; Hinshaw et al. 2003; Ko- 
matsu et al. 2003; Hinshaw et al. 2006; Spergel et al. 2006), 
testing the Gaussianity of the CMB has been imperative for 
our understanding of the Universe (Chiang et al. 2003; Gaz- 
tanagz & Wagg 2003; Coles et al. 2003; Park 2004; Eriksen 
et al. 2004b; Vielva et al. 2004; Cabella et al. 2004; Hansen 
et al. 2004; Mukherjee & Wang 2004; Larson & Wandelt 
2004; Naselsky et al. 2004; Tojeiro et al. 2006; Dineen & 
Coles 2003; Tegmark, de Oliveira-Costa & Hamilton 2003; 
de Oliveira-Costa et al. 2004; Eriksen et al. 2004a; Copi, 
Huterer & Starkman 2004; Schwarz et al. 2004; Land & 
Magueijo 2005; Bernui et al. 2006; Abramo et al. 2006; 
Chiang, Naselsky & Coles 2007; Cruz et al. 2006; McEwen 



in alphabetic order 



et al. 2006; Copi, Huterer Schwarz & Starkman 2006; Chi- 
ang, Coles & Naselsky 2007; Eriksen et al. 2007). 

One of the most general ways to test Gaussianity is 
based on the "random phase hypothesis" , as any departure 
from Gaussianity in the data shall register as some sort of 
phase correlation in the harmonic domain. There have been 
several non-Gaussianity methods devised from phase infor- 
mation: Shannon entropy of phases (Chiang & Coles 2000), 
phase mapping (Chiang et al. 2003), trigonometric moments 
(Naselsky, Doroshkevich & Verkhodanov 2004), phase sums 
(Matsubara 2003; Hikage et al. 2005), random walks (Stan- 
nard & Coles 2005; Naselsky et al. 2004), some of which 
have been deployed on WMAP full-sky maps and detection 
of non-Gaussianity has been decleared. 

As phases and morphology are closely related (Chi- 
ang 2001), one requirement for applying phases as a useful 
statistical diagnostic is continuity of the boundaries in the 
data, otherwise the phases would faithfully reflect boundary 
discontinuity by strong coupling. Therefore, those above- 
mentioned methods using phase information (particularly 
for CMB studies) can be deployed only on data with a full- 
sky coverage. 

Due to excessive foreground contamination near the 
Galactic plane, the WMAP science team has adopted a spe- 
cific foreground removal strategy using the so-called temper- 
ature masks (Bennett et al. 2003b; Hinshaw et al. 2006), 
which divide the full sky into 12 regions. The largest. Re- 
gion 0, covers about 89% of the full sky, whereas the other 
11 regions are masked due to heavy foreground emissions 
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Figure 1. The WMAP masks with the polar coordinates {6, if). 
The dark gray region is the WMAP Region 0, covering 89% of 
the full sky. The light gray area covers the WMAP Region 1 to 11 
(the Galaxy mask). The white dotted lines denote the Galactic 
latitude b = 28.18° (9 = 1.079) and -21.30° (6» = 1.943), the 
boundary of the Galaxy mask. 



of different kinds around the Galactic plane: synchrotron, 
free- free and dust emission (see Fig.l). Although a full-sky 
derived CMB map, the Internal Linear Combination (ILC) 
map, is combined from the 12 foreground-reduced regions 
and available to the public, most scientific results including 
the angular power spectrum are derived from the cleanest 
Region (Hinshaw et al. 2006), and the full-sky ILC map 
is known to still have foreground residuals near the Galactic 
plane. 

In this letter we present a new method for phase analysis 
on maps with Galaxy cut, assuming that the orthogonality 
of the Fourier series in the azimuthal direction outside the 
Galaxy cut is still preserved^. It is based on Fourier phases 
of equal-latitude pixel rings of the map, which is closely re- 
lated to the mean angle of the trigonometric moments on 
the full-sky phases with some weighting coefficients (Nasel- 
sky et al. 2004). We can examine the Fourier phases of all 
equal-latitude pixel rings from regions, e.g. WMAP Region 
0, while avoiding the polluted Galactic plane area. More im- 
portantly, we can pin down non-Gaussian features by using 
the phases derived this way, an advantage that is generally 
lacking in the analysis processed in harmonic domain. Note 
that all the above mentioned methods based on phases can 
be applied using the phases we derive in this letter. 



2 PHASES FROM AN INCOMPLETE SKY 
MAP 

The standard treatment for a full-sky CMB signal T{d, Lp) 
is via spherical harmonic decomposition: 



e=o m=-e 



(1) 



where i'max is the maximum multipole number used in map, 
9 and <p are the polar and azimuthal angle, respectively. 



^ Note that WMAP Region is not symmetric with respect to 
6 = 0, but |6| > 30° is surely outside the Galaxy mask (see Fig.l). 



and acm are the spherical harmonic coefficients. Y(m are the 
spherical harmonics, defined in terms of Legendre Polyno- 
mials: 



where 



NemPr'i^osd) exp{im(p), 



{21 + - my. 
47r(£-|-m)! 



(2) 



(3) 



The aim coefficients can be further written as acrn = 
\aern. \ exp[i$£m], where are the phases. If the CMB tem- 
perature anisotropies constitute a GRF, the real and imagi- 
nary part of the aim are both Gaussian distributed, or equiv- 
alently, the \aem\ are Rayleigh distributed and phases 
are uniformly random in [0, 27r]. In polar coordinate sys- 
tem 9 — tt/2 is associated with the Galactic plane {h = 0), 
as used by Healpix (Gorski, Hivon & Wandelt 1999) and 
Glesp (Doroshkevich et al. 2003) software packages. 

For signal from an incomplete sky coverage, implemen- 
tation of the spherical harmonic decomposition is no longer 
correct, as the orthogonality of the spherical harmonics Yem 
is broken (Gorski 1994). This is particularly the case when 
one is to analyze the WMAP ILC Galaxy-cut map. Nev- 
ertheless, Galaxy cut only breaks the orthogonality of the 
spherical harmonics over 9 direction, but not (p outside 
Galaxy cut (Gorski 1994). 

To see how phases of an incomplete sky map (e.g. ILC 
Galaxy-cut map) can be related to its full-sky phases, let 
us extract an equal- latitude pixel ring at ^ = ^c, where 
9c is outside the maximum latitude of any Galaxy masks. 
This ring T{9c,ip) = Tc(tp) is now one-dimensional signal, 
for which we can use a Fourier Transform approach with 
coefficients g^' 



Tc{<fi) = ^ 9m exp(imv5). 



where mmax ^ £max and 



3m = ^ / d<pTc{ip) exp{-iniip). 



(4) 



(5) 



We can then relate the ring with the full-sky signal via 
Eq.(l) and (2) and get 



ffm = 2^ NimPr {cos 9c) aer, 



(6) 



That is, the Fourier coefficients g^ of the ring can be 
expressed as a combination of the full-sky acrn- Writing 
ffm = |Sm|exp(i«;^), the phases are 



.1 E^|iV!i::i(^!l!^ 

Eotm! Wem{9c) cos^e 



(7) 



where Wem{9c) = NemPP {cos 9c)\aern\ - Note that the phases 
Km correspond to the "mean angle" of all with some 
weighting coefficients Wem{9c) involving the \aem\ (Nasel- 
sky et al. 2004). If the ring T{9c, (j)) is taken from a GRF, 
its phases combination of the uniformly random 

phases ^em, hence are also uniformly random in [0, 2tt]. We 
can then examine all the pixel rings of the ILC map for 
0^9^ 7r/3 and 27r/3 ^ 9 ^ tt without tapping the heavily 
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Figure 2. The mean angle (defined in Eq.(lO) and (11) with 
Am = 1 up to M = 50) of tlie Fourier phases from equal- 
latitude pixel rings Tc(ip) of the WMAP ILC 3-year map (top) 
and of a Gaussian realization (bottom). The gray area denotes 
the Galatic latitude boundary of the WMAP Galaxy mask at 
[—21.30°, 28.18°] (see Fig.l). One can see that the mean angles 
of the ILC map are fairly non-random, compared with the Gaus- 
sian realization. 



Figure 3. The mean angle (defined in Eq. (10) and (11) with 
Am = 1 up to M = 300) of the Fourier phases from equal- 
latitude pixel rings Tc((p) of the WMAP ILC 3-yeax map (top) 
and of a Gaussian realization (bottom). The gray area denotes 
the Galactic latitude boundary of the WMAP Galaxy mask at 
[—21.30°, 28.18°] (see Fig.l). One can see that the mean angles 
of the ILC map are significantly non-random, compared with the 
Gaussian realization. 



polluted region near the Galactic plane. Our demonstration 
here is a special case for a well known theory: any A'^ — n di- 
mensional cross sections of N dimensional Gaussian random 
process produce a Gaussian process as well. Thus, if one is to 
investigate the phases of the aem coefficients from a full-sky 
map, one can test alternatively the phases of equal-latitude 
pixel rings of the Galactic-cut map. 

However, a more intriguing question is whether we can 
reconstruct the phases of a full-sky signal ^irn by using the 
phases Km from the stripes of an incomplete sky map? Ob- 
viously we cannot reconstruct all the phases due to Galaxy 
cut, but we can recover significant paxt of the full-sky phases. 
Based on Gorski (1994) method and talcing into account 
that Galaxy cut map only breaks the orthogonality of the 
Legendre polynomials in 6 direction, there shall exist some 
polynomials KJ^(6) which are orthogonal to the Legendre 
polynomials Pi^{6) within some intervals [0, 7r/2 — ^cut] and 
[7r/2-|-6lcut,7r]. Namely, 

/ dxPr{x)K^{x) = Fie,m)Su', (8) 

*' ^CUt =COS &cut 



where F{i,m.) is the normalization coefficient. Then, defin- 
ing new coefficients 

SL = I dxgrJx)Kr{x) (9) 
= Nem.F{i,m)\aem\exp{i$ern); 

/ — ^cnt 
dxgm{x)KT{x) 

= NemF{i,m)\aem\exp(i^ern), 

which we can use for analysis of their phases. Since F{£, m) 
is a sign-flipping function, the phases of 5^ are equivalent 
to $^rn ± TT. Howevcr, the cross correlation of phases can 
be preserved. Care has to be taken in dcconvolution for the 
phases. Due to pixelization of the signal, particularly for the 
polar caps, modes at high multipole numbers tap the window 
function of the pixels. Implementing simple dcconvolution 
of the signal by window functions produces artifacts, which 
needs to be corrected by Tikhonov regularization. The same 
correction is needed for the high m modes as they are close 
to the Nyquist frequency. We will describe this approach in 
another paper. 
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3 MEAN ANGLE OF THE PHASES FROM 
THE ILC (GALAXY-CUT) MAP 

In this section, serving as an example of the Fourier phases 
K-m providing a useful diagnostic, we employ the trigono- 
metric moments and the mean angles on the phases derived 
from the equal-latitude pixel rings. The trigonometric mo- 
ments are defined as follows (Naselsky et al. 2004): 

M 

Cc(Am) = y^COs(w^+Am-Wm); 

m,= l 
M 

Sc(Am) = ^ sin (/t^+Am - Km) , (10) 

m— 1 

where M ^ imax — Am. Note that in this definition we use 
phase differences where Am ^ 1. The mean angle is defined 

as 

The mean angle can be seen as the resultant angle of 
Pearson's random walk (walks with a fix length in each 
step) : cxp[i(K^_,_Am ~k5^)] (Pearson 1906; Naselsky et 
al. 2004). For a GRF, the phases ^trn are uniformly random, 
so are the Km for each pixel ring. As the difference of any 
two random variables should be random as well, one then 
expects the mean angles G from an ensemble of Gaussian 
processes to be uniformly random in [0. 'liv]. 

We use the WMAP ILC 3-year map with /max = 512 
as an example of a high-resolution map. For each equal- 
latitude pixel ring Tc{ip), we use Fast Fourier Transform and 
obtain the phases k^- In Fig. 2 and 3 we plot the mean an- 
gles of each pixel rings with Am = 1 up to M = 50 and 
300 respectively against the Galactic latitude b. In each fig- 
ure, the top panel is the mean angles for ILC pixel rings 
and bottom for those from a Gaussian realization with the 
WMAP best-fit ACDM power spectrum. The gray area de- 
notes the Galactic latitude boundary of the WMAP Galaxy 
mask at [-21.30° , 28.18°] (see Fig.l). From both Fig.2 and 
3, one can see the ILC map outside the Galgixy mask has 
significant non-uniform distribution for the mean angles Q 
whereas for the Gaussian realization G are fairly uniformly 
random. Note that this example is for illustration purpose 
only, and more thorough analysis will be present in another 
paper. 



4 CONCLUSION 

In this Letter we have presented a new method of phase 
analysis of the CMB from an incomplete sky coverage. It 
is beised on Fourier phases of equal-latitude pixel rings of 
the underlying map, which are, theoretically speaking, re- 
lated to the mean angles of full-sky phases via well-defined 
weighting coefficients. We have also employed trigonometric 
moments and mean angles on the new phases, which has 
shown qualitatively significant non-random distribution of 
the mean angles, signature of departure of Gaussianity. We 
would like to emphasize that all the methods developed on 
using the full-sky phases can be easily implemented on the 
phases from an incomplete sky coverage. We will examine 



in details of non-Gaussianity using these phases in the next 
paper. 
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